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Abstract. 

Within the 1/N expansion of O(N) nonlinear a models for d < 4 it is possible to separate consis- 
tently the spin- wave and the massive-mode contributions to the scaling part of the free energy near 
criticality, and to evaluate them to 0(1 /N). For critical dimensions d = 2 + 2/n the Abe-Hikami 
anomaly is recovered, while for d = 2 the removal of the spin-wave term is justified. 



1 Hacettepe University, Physics Department, 06532, Beytepe, Ankara, Turkey 
fax: + (90) 312 2352550; e-mail: dilaver@thepl.phys.hun.edu.tr 

1 



Scaling is a fundamental property of physical systems in the neighborhood of a second order phase 
transition. In the field theory approach to the study of critical phenomena it is easy to identify 
the scaling properties of correlation functions. However it may be difficult to study scaling of 
bulk thermodynamical properties, and in particular the free energy. The reason behind this is 
the phenomenon of mixing: in a renormalizable field theory the vacuum expectation values of 
composite operators require in general not only a multiplicative renormalization, but also the 
subtraction of the contributions coming from lower-dimensionality operators carrying the same 
quantum numbers. These contributions, when present, are numerically dominant with respect 
to the scaling part of the expectation value, for the very simple reason that they are associated 
with lower powers of the (vanishing) mass scale, or inverse correlation length. Therefore one must 
in principle compute these terms with infinite accuracy (summing infinite orders of perturbation 
theory) in order for the evaluation of scaling contributions to become possible. In practice this 
" subtraction of perturbative tails" has been attempted with some success even in purely numerical 
computations of topological susceptibilities (l) . The reason behind this possibility stays in the fact 
that the physical degrees of fredom generating these contributions, the "spin waves", are quite 
different in nature from the (massive, topological) degrees of freedom associated with the scaling 
terms. Therefore, in a numerical simulation based on heating, the spin waves can be excited quite 
independently and earlier than the massive fields. When a plateau is temporarily reached, the 
value of the plateau can be nonperturbatively identified with the spin wave contribution, to be 
subtracted from the final (fully thermalized) vacuum expectation value. 

There are obvious conceptual and technical limits in the accuracy of this determination, and it 
would be nice to possess some analytical nonperturbative scheme in order to deal with this problem. 
In the context of two-dimensional O(N) lattice spin models, and of their regularized continuum 
counterparts, it was noticed some time ago that the 1/N expansion offered the possibility of 
consistently defining the "perturbative tails" of the free and internal energies, at least up to sec- 
ond nontrivial order ^). The procedure adopted might however appear to be tailored upon the 
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two-dimensional case, where criticality occurs in the zero-coupling limit, the models are pertur- 
batively renormalizable and asymptotically free, and the whole procedure essentially amounted 
to a prescription (principal value integration) for the evaluation of the integral representing the 
resummation of an infinite series of perturbative Feynman diagrams. We therefore decided to con- 
sider the problem of evaluating the scaling part of the free energy in O(N) nonlinear sigma models 
around criticality in dimensions 2 < d < 4, where criticality occurs at a finite value of the lattice 
coupling, and the critical exponents (hence the scaling properties) are nontrivial. 
Finding a subtraction procedure allowing for a nonambiguous identification of the scaling contri- 
butions cannot in this case be considered simply as the problem of making sense out of the formal 
sum of a known perturbative series, since we must first face the problem of correctly identifying 
which contributions come from spin wave degrees of freedom. 
For definitcness, let's consider the continuum version of nonlinear sigma models: 



S= ^N(3 J d d xd^.d^ (1) 
We label the coefficents of the 1 /N expansion for any physical quantity Q according to the notation: 



q = j2q* n ~* ( 2 ) 

i 

The (unsubtracted) free energy of the above model can be computed in the 1 /N expansion and its 
formal expression is ; 

N f d d n N 1 f d d n 

F = Y 7 (2^ ln/3(p2 + m ° ) ~T /3m o + 2 J J^- d H^-\ P ^)] + 0{l/N) (3) 
where mo is defined by the gap equation: 



d d p 1 
{2n) d p 2 + m 2 



P= I ,0-^-2 ■ ^2 (4) 



see ref H for cxplainations and notation. 
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We introduced the inverse propagator of the effective field: 



A-^mo) = \ j 



1 f d d p 



(2n) d q 2 + m 2 (p + q) 2 + m 2 
r(2-d/2) p 2 2 j/2 _ 2 d 1.3 p 2 

- 2 ( 47 r)^ (T +TOo) 2 ^ l(2 ~2'2'2'^T^ ) (5) 

where 2 F X (2 - |, |; |; p2 + m g ) is hypergeometric function f|. 

ft will be essential to our analysis that the inverse propagator admits an asymptotic expansion for 
small toq J5|, which is easily extracted from the representation: 



. , r(2 - d/2) r(d/2 - 1) 2 , 2 „ 2 x d - 3 (P-Pc) n/d-l , d 1 , 

A (p, TO ) = ^_ 2 ^ (p^^T+Jg-J^ 2 F 1 (_ J 1;_;___) ( 6 ) 

where (5 C is the (scheme-dependent) critical value of the coupling and one can show that, for small 

TOO, 

(^-^) = E |^K) d/2 - 1 
Such an asymptotic expansion can in turn be reinterpreted as the effect of an operator expansion. 
In order to identify the spin-wave contribution to the free energy we must evaluate the small-mass 
limit of eq.(^): 



FSW = f / (Sf 111 ^ + \ I ^A^m o) + O(l/N) (8) 



d<1 P 1 (R *\m 1 f d "P 

(27F ln(/3p) + 

where from eq.(|6|) we extracted the definition 



r(2-d/2)r(d/2-l) 2 r(l - d/2) to^ 2 

^' m ° j 2(47r) d / 2 r(d-2) P + (4tt)<</ 2 P 2 W 

The definitions eq.(||) and eq.(^|) were given with the aim of isolating those contributions to the free 

energy that are originated from a mixing with the identity operator. In order to justify our choice 

we must however prove that the subtracted free energy possesses the correct scaling properties. 



This is achieved by proving finiteness and scheme-independence of the dimensionless ratio; 



f = ±(F - F^)m- d (10) 

where we have introduced the mass gap m 2 (inverse correlation length), whose large N limit is m§. 
It is easy to obtain the 0(1 /N) correction in the form: 



m? 



The ratio / can also be expanded in powers of 1/N. The calculation of /o is straight-forward. The 
result is scheme- independent, and it is most easily derived by the use of dimensional regularization: 

h = \ [ j (0iMl + ^)-/3^Ho d = ^A- 1 (O,moK- d (12) 
The value of /i can be obtained from: 

h = \ [ S (0j ln [ A ob,m o )A- 1 ( P ,mo)]mo d - p ^ (13) 

In order to find a convenient representation for /i , we make the observation that the total derivative 
with respect to m\ can be represented in terms of the partial derivatives according to: 

d 9 d{3 d d , i . . d 

^ = ^ + 2A ~ (0 ' mo) ^ (14) 

Moreover the following identity holds under symmetric integration: 



1 d _■, 1 d— Id 1 

~^T^2 ^7^ A ' (P> m o)= „ 2 ^2 2F 1 ( — ,1;-; T —^-) (15) 



(p + im ) 2 + m 2 dp ^' uy p 2 + 4m2 1V 2 ' '2'i + J^ 

4m 



implying that: 

d d p 



m? = iA(0,m o ) J ^A(p,m )^A l {p,m ) (16) 
and as a consequence: 
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f l= l - m - d {J A^[lnA- 1 (p,mo)-lnA- 1 (p,m )-m^lnA- 1 (p,m )]} (17) 

We recognize that the expression appearing in the integrand of cq.(|l7]) must be O(mg), hence 
by naive power counting, confirmed by detailed analysis, the large-momentum behavior of the 
integrand turns out to be 0(l/p 4 ). This confirms our statement that fi is finite for all d < 4, 
independent of the regularization scheme adopted. 

Let's notice that / was evaluated for d = 2 in ref Q; and the result was: /o = ^ , /l = — 2/o ■ 
For d — 3 we obtained: fn = tti— 

J u 247T 



A 6 f°° 2 , n / 2 a\ , , 4 , 1 ,! 2 a; ,. , , 

— = — / x oi In (—arctan— )— In 1 — i—arctan - -) (18) 

fo TT J 1 V 2 7 1 7rx arctan^ v 2 a; 4 + x 2 ;J v ; 

S 1.97863 

Since for all d > 1 Aq" 1 vanishes for some real positive value of p 2 , we stress that our integration 
procedure must always be specified by the request that no imaginary part should appear in the 
final result (principal value prescription ) ](|. Our results for the range 1 < d < 4 are given by 
Figure 1. 

In order for our analysis to be complete, we must remember that in the critical regime even 
the high-frequency modes (spin waves) can bring a scaling contribution to the free energy Q. 
While logically distinct from the quantity we have just computed in the context of quantum field 
theories where a normal-ordering procedure can be defined in order to remove such a contribution, 
the scaling free energy originated from spin waves cannot be isolated in statistical systems and 
numerical simulations. Hence it is important to be able to compute the term proportional to m d 
in eq.(||). This is actually a reasonably simple task, since we can now exploit the properties of 
dimensional regularization and evaluate analytically eq.(^) for arbitrary d < 1 
From the explicit representation of Aq 1 ^, mo), eq.(||), we obtain in dimensional regularization: 
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1 f dd P ,„r^-2 , 4 r(d-2) , 2 . 



^ ^^[^j^^l-ra-^^) (19) 

where we have introduced the 0(1/ N) coefficient ai in the expansion of the critical exponent of 
the specific heat: 



d-4 2T(d) , on , 

00 = J— 2 ai = T(2 d/2)r(d/2)3 (20) 



It is important to notice that the imaginary part of the r.h.s. of eq.(19) is exactly needed in order 
to cancel the imaginary part originated by the naive integration of eq.(|l7|), which further justifies 
the procedure of removing all imaginary contributions. 

Here we face the phenomenon known as Abe-Hikami anomaly S; whenever 2/(d— 2) becomes a 
positive integer n, the expression in eq.(|l9|) is singular. However because of eq.(Q), in that same 
circumstance the expression we obtained is proportional to the positive integer power n+1 of (3—j3 c . 
We can therefore subtract a term analytic in (3 from eq. (|l^) and obtain the scaling contribution of 
spin waves to the free energy in the form: 

(sw) Nf m d 4 T(rf-2) » 

F ^ b ~2(d^T) [ d-2r(d/2~ir ] +0{1) (21) 

The most interesting feature of this result is the dependence on N, implying that for the special 
values d = 2+2 /n even the large N limit of the free energy is affected by the spin- wave contribution. 
When d — > 2, corresponding to in oo, singularities of the r.h.s. of eq.(20) accumulate on both 
sides of the real axis in the complex d plane. However it is easy to check that the limit exists when 
a small imaginary part is added to d, and its value is consistently 0, which further justifies the 
procedure of removing all spin-wave contributions in the two-dimensional case. 
The support from TUBITAK is acknowledged. 
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FIGURE CAPTIONS 



Figure 1. /1//0 versus d. 
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Figure 1: 
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